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ABSTRACT 


A technique has been developed for calculating feed-back and feed-forward 
gain matrices that stabilize a VTOL aircraft while enabling it to track input 
commands of forward and vertical velocity. Leverrier’s algorithm is used in 
a procedure for determining a set of state-variable, feedback gains that 
force the closed-loop poles and zeroes of one pilot-input transfer function 
to be at pre-selected positions in the s-plane. This set of feedback gains 
is then used to calculate the feedback and feedforward gains for the velocity- 
command controller. The method is computationally attractive since the gains 
are determined by solving systems of linear, simultaneous equations. Responses 
obtained using a digital simulation of the longitudinal dynamics of the GH-47 
helicopter are presented. 
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1. INTRODUCTION 


1.1 Background 

In recent years the overcrowding of this country's major airports has 
led to Increased interest in the development of commerical VTOL or STOL air- 
craft capable of operating in the 0-500 mile range. These vehicles would 
operate from separate runways at existing airports or from rooftops or short 
runways in or near business districts. Studies have shown (1,2) that such 
aircraft could reduce the total trip time providing they don't have to 
operate under existing take-off and landing procedures with their long 
delays. This will require innovations in vehicle design plus improved 
navigation, guidance, and control systems. One approach to this latter 
problem is to use on-board digital computers to handle the navigation, 
guidance, and control functions in an adaptive mode. Thus it is necessary 
to have a controller capable of following guidance commands from the computer 
and at the same time present acceptable flying qualities to pilot inputs. 

The Flight Instrxunentation Division of NASA's Langley Research Center 
is developing a technology base for digital adaptive control systems under 
the VTOL Advanced Landing Technology (VALT) program. This work was supported 
under this program by grant number NCR 47-018-005. 

1,2 Summary of Results 

This report presents the results of an investigation to develop a control 
algorithm as described above. The approach used in this investigation was to 
employ state variable feedback to achieve pole— zero placement of the vehicle's 
transfer functions. The initial attempt was to use the modified Newton- 
Ralphson technique to determine the necessary gains * This presented several 
difficulties, which are discussed later, and was abandoned in favor of a 



procedure, based on Leverrier's algorithm, that only requires the solution 
of linear simultaneous equa.tions. 

Using this procedure one first calculates a set of state-variable feed- 
back gains that result in acceptable responses to pilot inputs. This set 
of gains is then used to generate two sets of gains , one feedback and one 
feedforward, that maintain the same response to pilot inputs but enable the 
vehicle to track velocity commands from the guidance system. The algorithm 
is then used to control longitudinal dynamics in a digital simulation of the 
CH-47 helicopter during a decelerating approach. With this functioning 
properly, studies are conducted to demonstrate the controller's ability to 
function in the presence of errors in the estimates of aircraft parameters 
and at different sampling rates. 

1.3 Report Outline 

The equations and algorithms needed to express the coefficients of the 
closed-loop transfer function numerator and denominator polynomials are 
developed in Section 2 . In Section 3 the various methods used to calculate 
the feedback g ains are discussed along with some considerations in selecting 
a model. Section 4 presents the velocity-command configuration and procedure 
for calculating the necessary gains. In Section 5 the results of several 
fixed-point simulations are presented along with the simulation of a 
decelerating approach. Section 6 contains the results of the studies on 
varying aircraft parameters and sampling time. 
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2 . ANALYTICAL DEVELOPMENT 


2.1 Basic Concepts 

In considering the longitudinal dynamics of the CH-47 helicopter, the 


linearized equations of motion may be written in the form 
X = Ax + 


(2-1) 


where 

u (forward velocity perturbation) 
w (vertical velocity perturbation) 
q (pitch rate) 

6 (pitch angle perturbation) 


6 = 


6 (differential collective perturbation) 
e 

6 (collective perturbation) 
c 


A is the 4 x 4 differential transition matrix 
G is the 4x2 input matrix 

Introducing state -variable feedback from all states to both inputs 


results in a control vector of the form 


S = Kx + 6 

- - -p 


(2-2) 


where 

K is a 2 X 4 matrix of feedback gains 

6 is a 2 vector of pilot commands 
-p ■ ■ 

Applying this control to equation (2-1) yields 

X = (A + GK)x + G(S (2-3) 

— — “P 

The Laplace transform of equation (2-3) is 



This equation leads to eight closed-loop transfer functions relating the 
four state variables to the two controls. These transfer functions may be 
obtained by employing Cramer's rule. Thus, they all have the same denominator, 
I si - (A GK) I , but different numerators that are obtained by taking the 
determinants of the matrices that result when the appropriate columns of the 
G matrix are substituted for columns of [si - (A + GK) ] . Each transfer func- 
tion is a ratio of polynominals in s and the coefficients in these polynomials 
are functions of K. If values of K can be found to force these coefficients 
to take on pre-determined values, then the poles and zeroes of the closed- 
loop transfer functions can be placed anywhere in the s-plane. 

The determination of these coefficients requires the expansion of 4 x 4 
matrices where each element in the matrix is a function of the feedback gains . 
To do this by hand is an extremely tedious and error prone process. To avoid 
having to expand these determinants by hand, a procedure was developed, based 
on Leverrier’s algorithm, that calculates the desired coefficients. 

2.2 Leverrier's Algorithm 

Leverrier's algorithm is presented, with a proof, in reference (3) and 
used by Montgomery and Hatch (4) in their differential synthesis approach to 
pole-zero placement. The algorithm is an efficient computational technique 
for generating the coefficients of a system's open-loop transfer functions, 
or of the closed-loop transfer functions provided the feedback gains are 
known . 

In presenting the algorithm, consider equation (2-4) rewritten for an 
n X n open-loop system as 

x(s) = [si - (2-5) 

The term (si - A) ^ is of primary importance and may be written in the form 
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(si - A) 


-1 


^ + ... + B .s + B 
1 Z n-x n 


^ ^ ^ ^n-1® ^n 


n-2 


(2-6) 


where 


B . is an n X n matrix 

X 

d^ is a scalar coefficient of the characteristic equation. 

These B matrices and d coefficients are generated by Leverrier’s algorithm 
as presented below: 




®2 ■ V V 


= - trace (B^^A) 


= - 1/2 trace (B2A.) 




dj^ = - 1/k trace (B^^A) 


( 2 - 7 ) 


B = B ,A + d^ 
n n-1 n-1 


d = - 1 /n trace (B A) 
n n 


The transfer function between the state and input is of the 


form 


X. 

X 

6 . 

J 


ti-1 I ^ 

a. ..s + a. ,-s + ... + a _ 

ijl xj2 ijn 

n , . n— 1 I j n*— 2 , . j 

s + d.s + d„s + ... d ,s + d^ 

1 2 n-1 n 


(2-8) 


where the denominator is the same as in equation (2-6) and the a 's can be 
computed from 


a . = u. B, Gu. 

ijk —X k -3 


( 2 - 9 ) 


. T -I . -L .th j . th 

The u and u. are commensurable unit vectors wxth 1 xn the x and 3 posx- 

—1 -J 

tions respectively and 0 elsewhere. 

Thus the open-loop transfer function coefficients can be calculated, but 
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what is needed are the closed-loop coefficients expressed as functions of the 
feedback gains . 

2.3 Expansion of the Closed-Loop Denominator 

To achieve this consider the matrix [si - (A -I- GK) ] rewritten as 

Ksu^ " -1 ~ %^11 " -§2^21^ ■ -2 “ ~ -i-2^22^ 

, ( 2 - 10 ) 

(®H.3 “ £3 " £1^13 " £2^23^ ~ ^ ~ %^14 “ ^ 

where 

~i ^ column vector with 1 in the i position and 0 elsewhere 
th 

is the i column of the control matrix 
th 

a^ is the i column of the open-loop A matrix 

The determinant of this matrix is the closed— loop characteristic equation. 
Expanding this determinant using the rule 

l£i £2 ^-3 = 1^1 £2 ^3 ^4! + ^1 i£i £2 ££il <^2-ii) 

results in 


°CL ® 0 L 


2 

- I 

i=l j=l 


^ 4 4 X. ,6 

E 6^) + r ej 

X=1 J = 1 J c 


where 


( 2 - 12 ) 


= |aj^(s) a^Cs) ^^(s) a^(s) 


N(u,6^) = 

1&2. ^ ® ^ 

£3(3) a^Xs) 

u, 5 e 



,S > = 

1 % £2(3^ 

I2 —3 1 


is the open-loop denominator and H(x^,6j) 

open -loop transfer function. These terms are 

(5 

(2-7) and (2-9) . The N( ^ «^) , which will be 

*°c 

are not di re c t ly ob t ainab le us in g Leye rr ier ’ s 


is the numerator of the x./6. 

1 3 

easily generated using equations 
referred to as coupling functions, 
algorithm. 
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2 . 4 Determination, of the Coupling i^unctidns 


1 * e 

To develop the algorithm that generates N( ;; ) j consider the matrix 

1’ c 

formed by substituting the first column of the G matrix, times minus one, for 
the first column of the A matrix 

A' = [-g^ £2 £2 a^] (2-13) 

If Leverrier's algorithm is applied directly to the above A' matrix, one 


obtains 


d(s) 


[sV^ - A'] 


-1 




+ B, 
4 


(2-14) 

s 

^ + d^s^ + d2S^ + 

^ 3 ® + <14 


is 

an expression for 


-1 


'll 

“^21 "^13 

"^14 



■21 

^®“^22^ “^23 

“^24 


(2-15) 

bl 

"^32 ^®“^33^ 

-^34 



’41 

"%2 “^43 





where 



0 

0 

0 

0 


0 0 
1 0 
0 1 
0 0 


0 

0 

0 

1 


If B matrices and d coefficients, similar to those of equation (2-7), can 
he found for this matrix, then the coefficients of the coupling terms can be 
found by using equation (2-9) to solve a new system 

[sU^ - A']y^(s) - S .2 (2-16) 

tb 

where ^(s) is an n-vector of polynomials and all but the i element are 
coupling functions , 
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U^B^' -A'B2' = d^'I (2-19) 

U^B^' _ A-B3* = d3'I 
- A-B- = 

It was pointed out earlier that the d^s are merely the coefficients of N(u,6^), 
which are known, and can be determined. However, knowing d^’ and 
does not allow one to solve for B 2 * because is a singular matrix and there- 
fore has no inverse. 


Although the equations cannot be solved in the forward direction, they 

can be solved in reverse order if B. * can be determined and if is not 

4 

singular. In determining B^’ one must recall that it is the matrix of terms 
that do not contain s in the expression [sU^- A*] Now recall that B^ from 

equation (2-14) is the matrix of terms that do not contain s in the expression 

-I- ' ' 

[si - A’] . The only difference between [sU^- A’ ] and [si - A’ ] is that 
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one of the diagonal elements of [si - A'] is (s + while for [sU^- A'] 

it is just j • It is easy to show that the addition of s does not affect 

the matrix and therefore ® s matrix generated by Leverrier's 
algorithm. Knowing B^’ equation (2-19) can be solved in reverse order 

B 3 ’ = A’"^[U^B^' - d 3 'I] 

B^' = A'"^[U3B3' - d^'I] (2-20) 

Bj^' = - d^’I] 

The procedure may be summarized as follows : 

1. Apply Leverrier’s algorithm to the original A matrix and determine 

D = s + d-s + d-s +,*.+d-S + d 
ULi 1 z n^l n 

N(u, 6 ^), N(w, 6 ^), N(q, 6 ^), N(0,6^) 

N(u, 6 ^), N(w, 6 ^), N(q, 6 ^), N( 6 , 6 ^) 

2. Replace one column of the A matrix with and apply Leverrier's 

u ,6 u,(S u ,6 

algorithm to find B, . Determine B's and N( „^) , N( .^) , N(. ~^) . 

4 w, 6 ^ "q, 6 ^ ^ ^ 0 , 6 ^ 

3. Repeat step 2 for columns 2, 3, and 4 of the A matrix. 

2.5 Closed-^Loop Numerators 

To determine the closed loop numeratoi^ coefficients as functions of the 
feedback gains, Cramer's rule can be used in conjunction with equation (2-11), 
to yield 

( 2 - 21 ) 

3=1 3 c 

j#i 

J =1 J j’ e 
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Thus once the coupling function coefficients are knovm, both the closed-loop 
numerator and denominator polynomial coefficients may be determined as func- 
tions of K. 

2.6 Chen's Algorithm 

Since the development of the above algorithm for generating the coupling 
functions, the author has become aware of an algorithm developed by R.T.N, 
Chen (5) that will also generate these terms and which is more efficient 
computationally. The entire algorithm will not be developed here but the 
method of generating the coupling functions will be presented without proof. 
To simplify the expressions the following notation will be used; 



N(x^,6^) 


"n(x^.6^)“ 


N(x2.6e> 

N(x3,6^) 

1 A 
±2 = 

N(x2,6^) 

N(x3,6^) 


[n(x^,6^) 




Then the vectors of closed-loop numerators (using negative feedback) may be 
expressed as 


4 " il D~ 

CL ^OL 




(2-23) 


where is the vector of feedback gains. 

Comparing equations (2-23) and (2-22) it is apparent that the elements 

1 T X 

of the 4 X 4 matrix, the coupling functions . Since 

OL 

the diagonal elements of this matrix are zero there are three coupling 


functions per closed-loop numerator as previously determined. If one 
considers just one element of this matrix, for example the (1,2) element, 
it may be equated to a coupling function as shown below: 


S’«<= ■ ”0L 


(2-24) 
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The expression on the right hand side of the equal mark is the difference in 
two 6th order polynomials divided by a fourth order polynomial and it must 
yield a 2nd order polynomial with zero remainder because of the previous 
definition of the coupling functions. 

Since all of the coefficients of the KHS polynomials of equation (2-24) 
are directly obtainable using Leverrier's algorithm, generating the coupling 
functions merely requires a subroutine to multiply and divide polynomials. 
Furthermore, since it is known that the remainder must be zero, it is not 
necessary to calculate all seven coefficients of the products N(x^,6^) X 
N(Xj ,6^) and only the first three in the division. 

This approach has proven to be considerably more efficient than the 
previous method which involved four matrix inversions and many matrix multi- 
plications . 
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3 . DETERMINATION OF GAINS AND MODEL SELECTION 


3.1 Determination of Feedback Gains 

With the closed-loop transfer function coefficients expressed as functions 
of the feedback gains, the next step is to solve for a set of feedback gains 
that forces these coefficients to match predetermined values . This problem 
was approached from several different directions. 

Initially only the poles of the system were placed by forming a weighted 
least squares cost function and using nonlinear minimization techniques to 
match the coefficients. While this does place the poles, it has the obvious 
disadvantage of offering no control over the zero locations. It also yields 
non'-unique solutions since eight gains are available to match just four 
coefficients. 

A second approach was to introduce a 2 x 2 cross-coupling matrix and 
then specify two numerator polynomials as well as the denominator. Again 
the gains were obtained by forming a weighted least squares cost function 
and minimizing using a modified Newton-Ralphson technique. This approach 
was never completely satisfactory because it invariably resulted in large 
values of gain or unacceptable responses for the transfer functions not 
being specified. This is most likely caused by the fact that u, w, and 0 
for the helicopter are not independent. If two of these values are 
specified then the third must take on whatever value is required to 
satisfy the other two. The problem then lies in the fact that it is 
very difficult to a priori specify two transfer functions that yield 
a satisfactory response for the third variable. The above methods are 
discussed in more detail in (6) and (7). 

The method that was finally selected matches the poles and zeroes of 
one transfer function using only the feedback gains . Using this approach 
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it is only necessary to solve two systems of linear simultaneous equations. 

This method is presented by noting that the numerator coefficients are linear 
fiinctions of elements of the gain matrix. If desired zero locations in the 
s -plane are specified, then a model numerator polynomial can be determined. 
Equating the coefficients of like powers of s in the model and equation (2-22) 
yields a system of three linear simultaneous equations that can be solved for 
three elements of the gain matrix. Depending on whether a 6^ or 6^ transfer 
function is selected, three of the four elements of the first or second row 
of the K matrix are uniquely defined. The fourth element may now be specified 
and these four known values of gain substituted into equation (2-12). When 
this is done, the denominator coefficients become linear functions of the 
remaining four elements of the gain matrix. 

If the desired pole locations are now specified and a model characteristic 
equation formed, equating the coefficients of like powers of s results in a 
second linear system of simultaneous equations that may be solved for the 
remaining values of gain. The fact that the gains can be determined by 
solving two linear systems of simultaneous equations is a very attractive 
feature of this method and is discussed in more detail in reference (8) . 

3«2 Selection of a Model 

While, in theory, the above technique is capable of generating feedback 
gains that can place the closed-loop poles and zeroes of one transfer function 
at arbitrary locations in the s -plane, in practice, considerable care must be 
exercised in specifying these values and in selecting the extra value of gain 
from the numerator equations. If these values are not carefully selected, 
either the gains will be too large or the other transfer functions of the 
system will have zero locations that result in unacceptable responses. 

Examples of this are presented in Section 5. 
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Initially the pole and zero locations were selected to match specifica- 
tions set forth in references (9) and (10). This involved using pole-zero 
cancellation in order to match a second-order model and specifying two 
numerators. This effort was reported on in reference (11) but as stated 
previously, the general approach of trying to match two numerators did not 
produce satisfactory results and was abandoned. 

The final method used was basically trial and error in an attempt to 
satisfy the specifications shown in figure (1) . Wiiile no step-by-step 
procedure for selecting the pole-zero locations or the gain element can be 
given, the following general observations were made, 

1. Moving the pole locations farther out into the left-hand-plane to 
improve the speed of response requires that the corresponding zero locations 
also be moved out if the steady -state response is to be maintained. This, 
in general, will increase the gain magnitudes. 

2. In specifying the extra value of gain one should be aware of the 
relative magnitude of the state that the gain will multiply and the maximum 
allowable input to the system. This will allow an initial guess that is 
reasonable . 

3. Once the open-loop numerator, denominator and coupling fimctlon 
coefficients have been determined, it is quite easy to iterate on the pole- 
zero locations and extra gain to determine a set of acceptable values. 

4. For the CH-47, it required little modification to find acceptable 
values for all flight conditions considered, once an acceptable set of y 
values were determined at one flight condition. 

5. Eight-half -plane zeroes , while not desirable, will not necessarily 
cause problems if they are far enough out in the KEF or if the magnitude of 
the transfer function's response is small. 
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Figure (1) Typical Acceptable Responses To Velocity Commands 


4 . CLOSED-LOOP VELOCITY-COMMA.ND CONTROLLER 
Once the elements of the gain matrix have been determined so that the 
step input responses of all transfer functions are acceptable, then a satis- 
factory pilot-command control system of the form shown in figure (2) has been 
generated . 



Figure (2) System Configuration for Pilot-Command Controller 

With the above system it is possible, using the steady state responses of 
the transfer functions, to calculate the control inputs required to produce 
command changes in forward and vertical velocity (u^ and . However, this 
would be an open-loop controller with the inherent disadvantages of open- 
loop systems. To overcome this , the system configuration of figure (3) was 



Figure (3) System Configuration for Velocity-Command Controller 
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where 


10 0 0 
0 1 0 0 _ 

is a 2 X 2 matrix of feedforward gains 
is a 2 X 4 matrix of feedback gains 
c is a 2 vector of command changes in u and w 
The system dynamics are still defined by equation (2-1) , but the input becomes 


§. - ~ ^ 3 — (^- 1 ) 
§_ = K^x - K 2 K 2 X + K^c (4-2) 

^ - K^K^lx + K 2 C (4-3) 

or 

^ = Kx + K 2 C_ (4-4) 

Substituting into equation (2-1) yields 

X = (A + GK)x + GK 2 C (4-5) 


Now let K be the set of feedback gains obtained by specifying the poles 
and zeroes of one of the pilot-command transfer functions. Then it is 
apparent that the velocity-command and pilot-command transfer functions 
will have the same poles. Employing Cramer's rule and equation (2-11) it 
can be shown that the numerators of the u/u , w/u , u/w^, and w/w„ transfer 
functions may be expressed as 


N(u,u^) 


(4-6) 

N(w,u^) 


(4-7) 

N(u,w^) 


(4-8) 

N(w,w^) 


(4-9) 
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Each of the above equations represents a third order polynomial in s with four 

coefficients . While it is not possible to control all of the coefficients , it 

is possible to control one coefficient from each equation. If the system is 

to track velocity commands, the steady -state response of u and w to a unit 

step input on u^ must be one and zero respectively and the steady— state 

response of u and w to a unit step input on w_ must be zero and one respectively. 

This requires that the constant term of the numerator of the u/u transfer 

c 

function equal the constant term of the denominator (which is known) and the 

constant term of the w/u^ transfer function equal zero to provide de— coupling. 

Placing these constraints on equations (4-6) and (4-7) yields two linear 

simultaneous equations to solve for k and k„ . Similar restrictions on 

^11 ^21 

the constant terms of the w^ transfer functions result in two equations that 

yield k_ and k_ . Once the values of K„ are known K_ can be determined 
^12 ^22 ^ ^ 

using 

K = [K^ - K 2 K 3 ] (4_10) 
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5 . SIMULATION EESULTS 


This design procedure has been applied to the linearized longitudinal 
dynamics of the CH-47 helicopter. Results will be presented at several dif- 
ferent flight conditions and finally the results of a simulated decelerating 
approach will be presented. The first example will be for nominal velocities 
of 150 kts. forward and 250 fpm vertical descent. This example will be used 
to demonstrate the technique for calculating K 2 in the velocity -command 
configuration. 

5.1 Fixed Point I 

The A and G matrices shown below are based on stability and control 
derivatives provided by NASA's Langley Research Center. 


A = 


-0.05191 

0.02731 

-0.00014 

0 


-0.03898 8.8944 

-0.57793 250.67 

0.01769 -1.3044 

0 1.0 


-32.176 

1.2324 

0 

0 


G = 


-0.14909 

0.01857 

0.31973 

0 


-1.2698 

-8.9842 

0.22782 

0 


The open-loop characteristic equation for this A matrix is 


= s^ + 1.934s - 3 


2 

.579s - 


0.221s +0.0117 


The sign changes in the coefficients indicate an unstable system and right- 
half -plane poles occur at 0.034 and 1.2. 

Specifying the poles and zeroes of the w/6^ transfer function as 
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Poles = -0.75, -0.8, -0.8 ± jO.4 
Zeroes = -1.0, -0.8 ± jO.4 
resulted in the following gain matrix: 


0.0667 -0.02 -23.75 

-0.0021 0.0034 28.08 


•5.17 

0.324 


when specified at —0.02. This feedback gain matrix yields denominator 

and numerator polynomials for pilot-command transfer functions as shown below 
with the steady -state response of each listed to the right. 


°CL 

= 

+ 3.15s 

^ + 3.88s 

2 _ 

H 2.2s + 

0.48 

Steady -St ate 

Response/inch of control 

N(6^,u) 


0.15s^ 

“ 

9.23s^ 

- 

16 . Is 

- 6.33 

-13.2 

fps 

N(6^,w) 

= 

0.02s^ 

“ 

0,6s^ 

- 

s 

- 0.51 

- 1.06 

fps 

N(6g,q) 


0.32s^ 

+ 

0.21s^ 

H- 

0.01s 


0 


N(6^,0) 




0.32s^ 

+ 

0.21s 

+ 0 .01 

0.02 

rad 

N(6^,u) 


-1.27s^ 

“ 

10. 4s^ 

- 

11.4s 

- 1.58 

- 3.29 

fps 

N(6^,w) 

= 

-8.98s^ 


23. 4s^ 

- 

21.6s 

- 7.19 

-15.0 

fps 

N(Sj,,q) 


0.228s^ 

+ 

0.013s^ 

“ 

0.006s 


0 


N(6^,0) 




0.228s^ 

+ 

0 .013s 

- 0.006 

- 0.013 

rad 


Thus, the system has been stabilized and the pilot- command responses, 
shown in figures (4) and C5) , look reasonable. Sign changes in the u/6 , 
and numerators indicate right— half —plane zeroes ; however, in 

u/6^ and w/6^ they are located at 64 and 34 respectively and do not seriously 
affect the responses. The 0/S^ right-half -plane zero is at 0.14 and does cause 
some phase reversal. However, 0 is not very responsive to 5^ commands and 
this may not be objectionable. 

To determine the feedforward matrix, equations (4-6) - (4-9) are used 
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to form the following 2x2 systems of equations. 





— — 



-6.33 

-1.58 


^2 

^11 


0.48 

-0.51 

-7.19 


^2 

^21 


0 


and 


-6.33 

-1.58 


^2 

^^12 


0 

-0.51 

-7.19 


^2 

22 


0.48 


The resulting feedback and feedforward matrices are 



-0.0105 

0.0034 


-0.003 

-0.065 


-23.75 

28.08 


-5.17 

0.324 



-0.0772 

0.0054 


0.017 

-0.068 


These gains, used in the system configuration of figure (3), yield the follow- 


ing velocity -command transfer function numerators and steady -state responses; 


Steady— State 
Response to 
Unit Step Command 


N(u,u ) = -0.018s^ + 0.66s^ + 1.17s + 0.48 1.0 

G ■■■ . 

N(w,u^) = -0.050s^ - 0.08s^ - 0.04s + 0 0.0 

N(q,u^) = -0.023s^ - 0.02s^ - 0.0008s 0.0 

N(0,u^) = - 0.023s^ - 0.02s - 0.0008 0.002 rad 

N(u,w^) = 0.089s^ + 0.553s^ + 0.503s + 0 0.0 

N(w,w^) = 0.611s^ + 1.58s^ + 1.45s + 0.48 1.0 

N(q,w^) = -O.Ols^ + 0.003s^ + 0.0006s 0.0 

N( 6 ,w ) = — O.Ols^ + 0.003s + 0.0006 0 .001 rad 

c 
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The responses to velocity commands of = [10,0] and c^ = [0.10] were plotted 
and are shown in figures (6) and (7) respectively. The velocities settle to 
within 5% of the desired steady-state changes in 5 seconds with no overshoot 
and the decoupling between forward and vertical velocity commands is good. 

The pitch angle to vertical velocity command response shows a phase reversal 
that reaches a peak of 0.9 degrees in 1.3 seconds, but this may not be objec- 
tionable, If these responses were unacceptable to either pilots or passengers, 
the gains could easily be recalculated using new pole-zero specifications. 

5.2 Fixed Point II 

The results of the previous example do not completely satisfy the speci- 
fications of figure (1) in that the vertical velocity rise time exceeds two 

seconds. To correct this, the poles and zeroes of the w/6 transfer function 

c 

were moved to the following locations 
Poles = -0.9 ± jO.l, -0.9 ± jO.45 
Zeroes — -1.0, -0.9 ± jO.l 

k^2 was specified at -0.03. The nominal velocities were set at 110 fps 
forward and 500 fpm vertical descent since these values correspond to desired 
values along one section of the nominal approach path shown in figure (8). 

The A and G matrices will not be listed for the remaining examples, but can 
be determined based on the data in reference (12). The K and matrices 


for this example 

are shown below 



0.075 

-0.030 

-7.126 

-5.765 

K = 






-0.046 

0.050 

12.412 

2.156 


-0.062 

0.039 



K, - 





2 

0.007 

-0.108 
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The corresponding plots of the velocity responses to pilot commands and 
velocity commands are shown in figures (9) and (10) respectively. These 
responses do satisfy the previously stated specifications and the decoupling 
between u and w in the velocity-command mode looks good. There is some 
direction reversal in the u/w^ response but because of the small peak value 
of this response it should not be objectionable. The pitch and pitch rate 
responses are not shown since they were acceptable and are not of primary 
concern in this study. 

5.3 Fixed Point III 

This example is included to emphasize the importance of achieving satis- 
factory responses to pilot inputs prior to determining the gains for the 
velocity -command loop. Before attempting to simulate a decelerating approach, 
the controller’s ability to fianction properly at various nominal points 
throughout the flight range was investigated. In this example the nominal 
velocities were 20 fps forward and 3.33 fps vertical descent. The model 
and k ^2 were the same as in Example II. The velocity responses are shown 
in figures (11) and (12). The responses to 6^ commands shox^ extreme coupling 
between u and w and a rather severe direction reversal on u/6^ (again because 
of the small peak value this may not present a problem). The corresponding 
velocity-command responses demonstrate the same undesirable level of coupling 
between u and w for u^ commands. The peak value of 10 fps will most likely 
cause problems when flown in an adaptive mode. To correct this the value 
of ki 2 was adjusted (since it had been noted in previous studies that this 
can have considerable influence on the coupling between u and w) from -O.OS 
to +0.01. This yielded the responses shown in figures (13) and (14). The 
coupling has been greatly reduced and although the u/w^ response indicates 
a RHP zero the peak value is again small. 


26 








28 





Using the model from Example II and ” Q.Ol satisfactory responses 
were obtained at all flight conditions and these values were therefore used 
in the simulation of a decelerating approach. 

5.4 Decelerating Approach 

With the algorithm generating acceptable responses throughout the desired 
range of flight conditions, the next step was to check the controller's ability 
to function in an adaptive mode during a decelerating approach. 

The section from C to D of the nominal approach path of figure (8) was 
used for this study with the following modifications. 

1. The vertical velocity between B and C was set at 0 rather than 500 fpm. 
(The controller was checked under both conditions but the results presented 

are for the latter since this is a more severe test of its ability to come onto 
and remain on the desired path.) 

2. The vertical velocity was linearly decreased from 500 fpm to 0 over 
the range from 100 feet to 0. This takes both forward and vertical velocities 
to zero simultaneously and brings the vehicle to hover at an altitude of 50 
feet. 

In order to simulate a decelerating approach a guidance law had to be 
developed that would command the changes in forward and vertical velocities 
necessary to keep the vehicle on the desired approach path. This requirement 
is complicated by the fact that the controller is designed to command changes 
in body axis velocities while the approach path is defined in terms of approach 
navigation frame (ANF) velocities. The relationship between these two axis 
systems is shown in figure (15) . 
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Body Axis u 



^ ANF X 


Body Axis w 

Figure (15) Relationship Between Body Axis and ANF Axis 


These systems are related by the equations 
u = X cos 0 - z sin 0 

(5-1) 

w = X sin 0 + z cos 0 

Thus the basic procedure is to determine the desired changes in x and 
z and then use equation (5.1) to generate the necessary commands in u and w. 
The equations for generating x and z are of the form proposed by Hoffman, 
Zvara, Bryson, and Ham (13). The general procedure is to develop expressions 

.» ft ' ' 

for X, z, and z as functions of range and determine the necessary purturba- 
tions with the following equations : 

Ax = - ^aCTUAL^ “ ACTUAL “ ^DESIRED^ 

(5-2) 

^ ^ ^DESIRED “ ^ACTUAL^ ^z ^ ^ACTUAL ~ ^DESIRED^ 

where c and c are weighting factors whose function is to maintain the 
X z 

desired altitude profile. For the results to be presented c^ = 0.5 and 
c^ =1.5. The equations expressing x, z, and z as functions of range are 
given below and were obtained by fitting curves to desired profiles of 


figure (8) . 


For 2500 < R <1200 


Desired 

^desired 

^DESIRED 


0.0333 • R + 26.0 
8.333 

-1.21212 • R^+ 0.1448484 ' R + 173.636 


(5-3) 
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For 1200 < R < 400 


Desired “ • r + 12.0 

^DESIRED “ ^*333 

Zjjesired “ ■*■ ^ 

For 400 £ R < 0 

Desired ~ ^ + 0.12 • r 


^DESIRED 


8.3333 

0.08333 


400 £ R £ 100 
100 < R < 0 


^DESIRED 


+ 0.55 • R + 50.0 


In the ^j)£gjg£j) equations the value of R is the actual value at the time 
the velocity commands are updated. For ^gglRED ^DESIRED value of R 
is obtained by subtracting DT • x^cTUAL’ where DT is the guidance update time, 
from the actual value of R. This does not give an exact value for R and 
therefore x and z at the next update point since x is changing but it does 
yield an estimate at some future point along the trajectory and it was 
reasoned that this is better than commanding x and z based on the current 
location. 


For the results to be presented the sampling rate was 10 per sec. and 
the guidance commands were updated every 0.5 sec. while the actual aircraft 
parameters and gains were updated every second. The model and value of 
were the same as the last example at fixed point III. 

The results of this simulation appear la figures (16) through (20). The 
plot of altitude vs. range shows good tracking while the ANF x vs. range curve 
shows the actual velocity consistently above the desired value. It should 
be emphasized, however, that the guidance law, equation (5.2) penalizes errors 
in altitude so this result might have been anticipated. It should also be 
emphasized that this guidance law was rather hastily devised to test the 
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50.5. 

45,5- 




Figure (16) ANF Forward Velocity va. Range 
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feasibility of using the controller jn an adaptive decelerating mode and 
it has demonstrated that this is possible. The remaining plots of MF z, 
pitch rate, and pitch angle vs. time are included to demonstrate that these 
values lie within acceptable limits. 

Thus the controller appears capable of performing the fimctions it was 
designed for without exceeding the capabilities of the aircraft. The item 
remaining to be investigated is the controller's ability to function in the 
presence of errors in the estimates of aircraft parameters and at different 
sampling rates. The results of this study are presented in the next section. 
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6, SENSITIVITY STUDIES 


Since the control algorithm is to be used in an on-board computer where 
the aircraft parameters must either be stored or estimated, it is necessary 
to have information concerning the controller's ability to function in the 
presence of errors in these parameters. The computer must also operate in 
real time and perform several functions in addition to the control function. 
Therefore it is necessary to know the controller’s sensitivity to varying 
sampling rates. The results of these studies appear in this section. 

6.1 Sensitivity to Parameter Errors 

The linearized equations that describe the vehicle’s motion about a 
nominal point are of the form of equation (2.1) which is written below in 
expanded form to show the actual entries in the A and G matrices. 


- - 


— 




— 


. 


— 

— 

u 



^w/m 


-g cos 

0 

n 


u 


e/m 

c/m 

w 

= 

Vm 

^w/m 


-g sin 

0 

n 


w 


e/m 

^6 / 
c/m 

q 


\/l 

yy 

\/l 

yy 

“q/l 

yy 

0 



q 


^e/I 

yy 

^c/l 

yy 

0 


0 

0 

1 

0 



0 


0 

0 


where 



= nominal forward velocity 
= nominal vertical velocity 
0^ > nominal pitch angle 
g == gravitational acceleration 

The remaining terms in the matrices are the stability and control derivatives 

that must be stored or estimated. The procedure used in the program was to 

specify U and W and use these values to determine 0 and the derivatives 
" n n n 

based on the data in reference (12) . 
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To determine the effect of errors in these parameters, values of and 
were specified and the necessary values of gain were calculated, using the 
techniques previously discussed. These gains, based on correct values of 
the vehicle’s parameters, were then used to control the vehicle with the 
parameters changed, one by one, by multiplying by 0 and 3. This yields a 
fairly wide range of variation and should indicate the relative importance 
of each parameter. 

The corresponding responses in forward and vertical velocity were plotted 
for velocity commands of u^ = 10 fps and w^ = 0 fps for nominal velocities 
of U == 110 fps and W 8,33 fps. The sampling rate was 10 per second. 

These responses appear in figures (21) through (36) where each plot contains 
three curves representing the response, with the correct value, the response 
with the indicated parameter set to zero and the response with the parameter 
multiplied by three. 

Based on these curves it appears that the controller is fairly insensitive 
to variations in the vehicle parameters. For the stability derivatives, 
variations in the Z and M derivatives show considerable influence over the 
decoupling between u and w in steady state, but very little d..fference in 
the response for the first one or two seconds. If the controller is being 
flown in an adaptive mode, this may be the critical portion of the response 
since commands in u and w will most likely be updated every 0.5 to 1.0 


seconds . The control derivatives 
become unstable when set to zero. 


Zp y and /I cause the vehicle to 
This however might have been expected 


since setting these parameters to zero severely limits the controller’s 


ability to influence the vehicle’s motion. The control derivative 


Mjs /I , when multiplied by three, also caused the vehicle to become 
6^ yy’ 

unstable. Howevetj a factor of three error in a control derivative is 
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Figure (25) Velacity Responses to u = [10,0] With X 
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an extreme variation. 

The above results are admittedly limited in their application since the 
vehicle parameters were varied one by one. They should however serve to 
indicate the relative importance of the parameters. It is recommended that 
work in this area be continued. 

6.2 Sensitivity to Sampling Rate 

Using the procedures described in this report the gains are calculated 
assuming a continuous system with dynamics described by equation 2.1, 
rewritten below. 

x = Ax + G^ (2.1) 

Since the actual system, however, is digital, the equations of motion become 

^i+l = ^ + B ^ (6.1) 

where 

T 

F = / AF dt F = I (6.2) 

o 0 


T 

B = / (AB + G) dt B = 0 
^ o 


(6.3) 


with T being the sample time. Thus the procedure was to calculate the gains 
for a continuous system and use those gains to generate ^ in equation (6.1) 
using the following equation 


6 ... + 1 = K X. 
— x+1 — X 


+ K„ 6 , 

2 “C l 


(6.4) 


Equation (6.4) can represent the input for the velocity command controller by 

setting equal to the feedforward gain matrix and letting be a vector 

of velocity commands or it can represent the input for pilot commands by 

setting equal to an identity matrix and letting 6 , be the vector of 
/ — Cl 

pilot inputs . 
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Applying the input of equation (6.4) to equation (6.1) the state responses 

were plotted for different sampling rates by generating new F and B matrices 

using equations (6.2) and (6.3). Responses were plotted for values of T 

between Q.05 and 2.0 sec at nominal velocities of U =20 fps, W = 3.33 fps 

n ^ ’ n 

and U = 110 fps, W =8.33 fps. 

For nominal velocities of = 20 fps and = 3.33 fps there is very 
little change in the vehicle’s response for values of T of 0.5 second and 
smaller. The responses for T = 1.0 second were still stable but were 
becoming oscillatory and at T = 2.0 seconds the system was unstable. At 
the higher velocities of = 110 fps and = 8.33 fps the responses had 
become oscillatory at T = 0.5 second and unstable at T =1.0 second. For 
T = 0.33 second and below, however, there was little difference indicating 
that these rates should be acceptable at all points along the approach 
trajectory . 
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7 . CONCLUSIONS 


A technique has been developed for generating feedback and feedforward 
gains capable of producing acceptable vehicle responses to both pilot and 
guidance system commands. The calculations are simple and lend themselves 
to on-board computation in an adaptive mode of operation. Sensitivity 
studies indicate that the system can ftmction in the presence of errors in 
vehicle parameters and over a fairly wide range of sampling intervals. 

Publications 

The results of this work have been published as listed below in 
reference 8 and in the paper "A Velocity Conimand Controller For A VTOL 
Aircraft" presented at the 1976 JACC, West Lafayette, Indiana, July 1976. 
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A TECHNIQUE FOR POLE-ZERO PLACEMENT 
FOR DUAL-INPUT CONTROL SYSTEMS 


Gerald F. Reid 
Virginia Military Institute 
Lexington, Virginia 


ABSTRACT 

An algorithm is presented for determining the feedback gains required to 
place the poles and zeros, of pre-selected transfer functions of a dual-input 
control system, at arbitrary locations in the S-plane. The technique employs 
Leverrier's algorithm and has been used in a digital simulation to control 
the longitudinal dynamics of the CH-47 helicopter. 
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A VELOCITY -COMMAND CONTROLLER EOR A VTOL AIRCRAFT 


Gerald F. Reid 
Virginia Military Institute 
Lexington 5 Virginia 


ABSTRACT 

A technique is presented for calculating feedback and feedfon^rard gain 
matrices that enable a VTOL aircraft to track input commands of forward and 
vertical velocity while maintaining acceptable responses to pilot inputs. 
Leverrier’s algorithm is used for determining a set of state-variable, feed- 
back gains that force the closed-loop poles and zeros of one pilot-input 
transfer fimction to pre-selected positions in the s-plane. This set of 
feedback gains is then used to calculate the feedback and feedforward gains 
for the velocity-command controller. The method is computationally attractive 
since the gains are determined by solving systems of linear, simultaneous 
equations. The method has been used in a digital simulation of the CH-47 
helicopter to control longitudinal dynamics . 
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